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ABSTRACT 

Using a suite of progenitor models, neutrino luminosities, and two-dimensional (2D) simulations, we inves- 
tigate the matter gravitational-wave (GW) emission from postbounce phases of neutrino-driven core-collapse 
supernovae (CCSNe). These phases include prompt and steady-state convection, the standing accretion shock 
instability (SASI), and asymmetric explosions. For the stages before explosion, we propose a model for the 
source of GW emission. Downdrafts of the postshock-convection/SASI region strike the protoneutron star "sur- 
face" with large speeds and are decelerated by buoyancy forces. We find that the GW amplitude is set by the 
magnitude of deceleration and, by extension, the downdraft's speed and the vigor of postshock-convective/S ASI 
motions. However, the characteristic frequencies, which evolve from ~100 Hz after bounce to ~300-400 Hz, 
are practically independent of these speeds (and turnover timescales). Instead, they are set by the deceleration 
timescale, which is in turn set by the buoyancy frequency at the lower boundary of postshock convection. Con- 
sequently, the characteristic GW frequencies are dependent upon a combination of core structure attributes, 
specifically the dense-matter equation of state (EOS) and details that determine the gradients at the boundary, 
including the accretion-rate history, the EOS at subnuclear densities, and neutrino transport. During explosion, 
the high frequency signal wanes and is replaced by a strong low frequency, ~ 10s of Hz, signal that reveals the 
general morphology of the explosion (i.e. prolate, oblate, or spherical). However, current and near-future GW 
detectors are sensitive to GW power at frequencies >50 Hz. Therefore, the signature of explosion will be the 
abrupt reduction of detectable GW emission. 

Subject headings: hydrodynamics — instabilities — shock waves — supernovae: general — gravitational 
waves — dense matter — equation of state — turbulence 
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1. INTRODUCTION 

Core-collapse supernovae (CCSNe) are among the most en- 
ergetic events in the Universe; they herald the birth of neutron 
stars and black holes, are a major site for nucleosynthesis, in- 
fluence galactic dynamics, trigger further star formation, and 
are prodigious emitters of neutrinos and gravitational waves. 
Hence, it is important to understand the mechanism of explo- 
sion, yet the details have remained elusive for many decades. 

Gravitational waves (GWs) promise to be important di- 
agnostics for revealing the secrets of the CCSN mecha- 
nism. Despite the elusiveness of the CCSN problem, current 
multi-dimensional simulations show several promising mech- 
anisms. Whether the solution is the convection/SASI-aided 
neutrino mechanism (Marek & Janka 2009; Murphy & Bur- 
rows 2008b), the acoustic mechanism (Burrows et al. 2006, 
2007c), or magnetohydrodynamic (MHD) jets (Burrows et al. 
2007b; Dessart et al. 2008 and references therein), viable 
explosion mechanisms for all but the least massive of mas- 
sive stars (Kitaura et al. 2006; Burrows et al. 2007a) appear 
fundamentally multi-dimensional. Because photons couple 
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strongly with matter, the mechanism is obscured by many so- 
lar masses and traditional means of observation are limited. 
However, neutrinos and GWs interact very weakly with mat- 
ter and provide a clear view into the dynamics of the mecha- 
nism. Unlike neutrinos, GWs are emitted (at lowest order) by 
time-changing quadrupole mass/energy motions. Therefore, 
the multi-dimensional mechanisms for CCSNe are expected 
to be a significant source of GWs, which will in turn provide 
a diagnostic into their multi-dimensional nature. 

This important connection between CCSNe and GWs was 
recognized quite early (for a comprehensive review of GWs 
and CCSNe, see Ott 2009b). Most early investigations were 
limited to the collapse and bounce of rapidly rotating iron 
cores (Miiller 1982; Monchmeyeret al. 1991; Yamada & Sato 
1995; Zwerger & Miiller 1997; Dimmelmeier et al. 2002; 
Kotake et al. 2003; Ott et al. 2004). Current investigations, 
that focus on these early phases and rapid rotation, use 2D 
or 3D simulations, conformally-flat or full general relativ- 
ity, finite-temperature nuclear equations of state (EOS), and 
deleptonization during collapse (Dimmelmeier et al. 2007, 
2008; Ott et al. 2007). However, the rotation rates neces- 
sary for GWs with sizeable amplitudes are likely to obtain in 
only ~1% of the massive star population (Heger et al. 2005; 
Hirschi et al. 2004; Woosley & Heger 2006; Ott et al. 2006b). 

Another set of studies has focused on the GW signatures of 
postbounce phases and dynamics. These include convection 
in the protoneutron star (PNS) and in the postshock region 
(Burrows & Hayes 1996; Miiller & Janka 1997; Miiller et al. 
2004; Marek et al. 2009), the standing accretion shock insta- 
bility (SASI) (Kotake et al. 2007, 2009; Marek et al. 2009), 
PNS internal g-modes associated with an acoustic mechanism 
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for explosion (Burrows et al. 2006, 2007c; Ott et al. 2006a), 
rapid rotation in the context of MHD jets (Kotake et al. 2004; 
Obergaulinger et al. 2006), or asymmetric explosion and neu- 
trino emission which lead to a long-term nonzero GW strain, 
or "memory" (Braginskii & Thorne 1987; Burrows & Hayes 
1996; Miiller & Janka 1997; Miiller et al. 2004). Each sig- 
nature is associated with the three possible explosion mecha- 
nisms, and Ott (2009a) suggests that their GW characteristics 
may help to distinguish which mechanism obtains in a partic- 
ular supernova. 

In this paper, we do not address the GW signatures of PNS 
g-modes nor rapid rotation but focus on the GWs from the 
postbounce phases of the neutrino mechanism. The bounce 
shock is quickly stalled by nuclear dissociation, electron cap- 
ture, and neutrino losses (Mazurek 1982; Bruenn 1985, 1989) 
and leaves a negative entropy gradient. This leads to an initial, 
short-lived (< 50 ms) phase of strong convective overturn, 
prompt convection (Burrows & Fryxell 1992), and an asso- 
ciated burst of GWs (Marek et al. 2009; Ott 2009b). More 
than two decades ago, Wilson (1985) and Bethe & Wilson 
(1985) suggested that a fraction of the neutrinos being emitted 
from depth (<100 km) would be recaptured in the gain region 
(>100 km), reviving the stalled shock into explosion. How- 
ever, detailed one-dimensional (ID) simulations have shown 
that this mechanism, the neutrino mechanism, fails in ID 
(Liebendorfer et al. 2001b, a; Rampp & Janka 2002; Buras 
et al. 2003; Thompson et al. 2003; Liebendorfer et al. 2005), 
except for the least massive of the massive stars (Kitaura et al. 
2006; Burrows et al. 2007a). Recent 2D simulations that 
are subject to aspherical instabilities, specifically postshock 
convection, which is driven by neutrino heating in the gain 
region and the SASI, suggest that the neutrino mechanism 
may yet succeed, though it fails in ID (Herant et al. 1994; 
Janka & Miiller 1995; Burrows et al. 1995; Janka & Miiller 
1996; Burrows et al. 2007c; Kitaura et al. 2006; Buras et al. 
2006b,a; Marek & Janka 2009; Ott et al. 2008; Murphy & 
Burrows 2008b). In particular, Murphy & Burrows (2008b) 
parametrized the neutrino luminosity and found that the criti- 
cal neutrino luminosity for explosions is lower in 2D simula- 
tions compared to ID. These results suggest that asymmetries 
are not only important for the production of GWs, but also for 
the success of the neutrino mechanism. 

PNS convection (see Dessart et al. (2006) and the refer- 
ences therein) and postshock convection (Burrows et al. 1995; 
Janka & Miiller 1996) have been recognized as an important 
ingredient in CCSNe for many decades. However, it was only 
recently that the SASI has been recognized as a distinct in- 
stability that augments neutrino-driven convection in the gain 
region (Blondin et al. 2003). The mechanism responsible for 
this instability is debated to be either an advective-acoustic cy- 
cle (Foglizzo & Tagger 2000; Foglizzo 2001, 2002; Blondin 
et al. 2003; Ohnishi et al. 2006; Foglizzo et al. 2007; Iwakami 
et al. 2008b; Scheck et al. 2008; Yamasaki & Foglizzo 2008; 
Foglizzo 2009; Sato et al. 2009) or purely acoustic cycle 
(Blondin & Mezzacappa 2006; Laming 2007), though recent 
linear and nonlinear analysis of simple models suggest that 
the advective-acoustic cycle is responsible (Sato et al. 2009; 
Foglizzo 2009). An additional complication is that it is hard 
to disentangle the effects of the SASI from those of postshock 
convection. They both occupy the same region, reach nonlin- 
ear saturation quickly, and likely influence each other's non- 
linear motions. Given the recent recognition of the SASI, it 
is no wonder that few have investigated fully the GW char- 



acteristics of the convection/SASI-aided neutrino mechanism 
(Kotake et al. 2007, 2009; Marek et al. 2009). 

Using simulations performed with the 2D radiation- 
hydrodynamics code VULCAN/2D, Ott et al. (2006a) de- 
scribed some GW signatures of the SASI, but primarily fo- 
cused on the GW emission from strong core g-modes. Ko- 
take et al. (2007) were the first to specifically investigate the 
GW signature of the SASI in 2D. To expedite calculations, 
they made several assumptions. Their initial conditions were 
derived from a steady-state approximation of the postshock 
structure and not from stellar evolution calculations. They 
used an inner boundary at 50 km and kept the density at this 
boundary constant at 10 11 g cm -3 . Although they used a 
relativistic-mean-field EOS (Shen et al. 1998), they approx- 
imated the effects of neutrino interactions with local heating 
and cooling terms and assumed a fixed accretion rate at the 
grid outer boundary. Even though they omitted the PNS and 
proper neutrino transport, they estimated the GW emission 
from matter and neutrinos, and found that the GW amplitudes 
due to asymmetric neutrino emission are ^100 times that of 
matter, but at low frequency and hence lower power. They 
also found that the matter GW signal grows as the SASI en- 
ters the nonlinear phase and that the characteristic frequency 
from the matter GW emission is ~100 Hz. Interestingly, 
they did not find a signature of explosion that is distinct from 
steady-state convection and/or SASI. More recently, Kotake 
et al. (2009) used a similar simulation approach in 3D simu- 
lations and employed a ray-tracing technique to estimate the 
asymmetric neutrino emission, obtaining a neutrino GW sig- 
nal comparable in magnitude to the matter signal. Though 
one might expect higher neutrino luminosities to result in 
stronger convection/S ASI motion, they found that the strength 
of the GW emission does not correlate with neutrino luminos- 
ity. Rather, they concluded that the stochastic motions of the 
SASI prevents such a correlation between GW strength and 
neutrino luminosity. 

Marek et al. (2009) use different approximations to in- 
vestigate the GW characteristics in 2D simulations and ob- 
tain different results. They use a compressible liquid-drop 
model EOS as well as relativistic-mean-field and a Brueckner- 
Hartree-Fock EOSs, ID ray-by-ray neutrino transport, ini- 
tial conditions derived from stellar evolution calculations, a 
pseudo-GR gravitational potential that mimics the effects of 
GR in spherical symmetry, and an inner boundary at ^2 km. 
With these more realistic assumptions, they obtain matter GW 
amplitudes that are only half of the neutrino GW signal. Just 
after bounce, they report characteristic frequencies of ~400 
Hz, similar to Kotake et al. (2007). However, Marek et al. 
(2009) associate these frequencies with the short-lived prompt 
convection rather than general SASI motions. When non- 
linear SASI begins in earnest, augmenting the nonradial mo- 
tions of neutrino-driven convection, they report that the GW 
power peaks at ^300-800 Hz and attribute the source of GWs 
to vigorous SASI and convective motions above and below 
the neutrinosphere. A primary conclusion of their work is 
that the peak frequency depends upon the compactness of the 
PNS and, by extension, the dense-matter EOS, in that softer 
EOSs result in higher frequencies. 

Both groups state that GWs are generated by convective 
overturns and sloshing motions of the SASI. In fact, Kotake 
et al. (2009) make direct reference to overturn and sloshing 
timescales for the origin of GW characteristic frequencies. 



3 



Marek et al. (2009) also discuss these timescales in the con- 
text of prompt and PNS convection, but they go further to 
suggest that asymmetric motions in the neutrino heating and 
cooling layers stirred by vigorous SASI funnels and convec- 
tive overturn are significant sources of GW emission as well. 
However, neither quantitatively define the relevant timescales, 
nor give support for these associations. In fact, most timescale 
definitions in the literature associated with postshock convec- 
tion and the SASI give frequencies (~30-100 Hz) that are 
quite low (Burrows et al. 2007b; Scheck et al. 2008; Marek & 
Janka 2009) compared to the frequencies (>300) associated 
with peak GW power. This discrepancy suggests that another 
timescale is relevant in determining the characteristic GW fre- 
quency. So, here we ask several obvious questions. What de- 
termines the characteristic frequencies and amplitudes, how 
do these change with progenitor mass and neutrino luminos- 
ity, and what is the GW signature of explosion? 

Theoretical GW signals that reproduce observations with 
fidelity will require 3D radiation-hydrodynamic simulations. 
However, such simulations are computationally expensive and 
prohibit systematic studies of the important physics. In the 
spirit of Murphy & Burrows (2008b), we simulate simplified 
models that nevertheless retain the important physics, and al- 
low one to adjust key parameters. Our approach is a compro- 
mise between Kotake et al. (2007) and Marek et al. (2009). 
Like Marek et al. (2009), we use recent core-collapse progen- 
itor models (Woosley & Heger 2007) as initial conditions and 
a finite-temperature relativistic-mean-field EOS. Rather than 
studying the effects of the EOS, we perform simulations with 
four progenitor models: 12, 15, 20, and 40 M Q . To avoid 
waiting six months (or more!) for a single simulation, we 
compromise on the neutrino transport and use local heating 
and cooling approximations similar to Kotake et al. (2007). 
Finally, we use spherical Newtonian gravity. Though these 
assumptions prevent accurate models, simulations can be cal- 
culated in a timely manner that include important qualitative 
features of core collapse. Because we do not use neutrino 
transport, but a simple local heating and cooling algorithm, 
we calculate the GW signal due to asymmetric mass motions 
and not asymmetric neutrino emission. Since the neutrino and 
matter GW signals are emitted at low and high frequencies, 
respectively, and GW detectors are sensitive to the high fre- 
quencies, calculating the matter GW signal is sufficient for 
predicting GW observations in the near future. 

In the following sections, we present a systematic 
parametrization of progenitor mass and neutrino luminosity 
in the neutrino mechanism context. In §2, we present the 
basic equations, assumptions, and numerical techniques. In 
§3, we present and discuss the results of this paper. For GW 
extraction, we consider the quadrupole formula of the slow- 
motion, weak-field approximation (§3.1). In §3.2, we de- 
scribe the generic features of GW emission, which include 
prompt convection, postshock convection, SASI, and asym- 
metric explosions. In §3.3, we show energy spectra and spec- 
trograms, describing the characteristic frequencies and their 
dependence on progenitor mass and time of explosion. In 
§3.4, we consider the characteristic frequencies and ampli- 
tudes of the downdrafts striking the PNS "surface," where 
they decelerate due to buoyancy force and show that these 
match the characteristic frequencies and amplitudes of the 
GWs. Finally, in §4, we summarize our results and suggest 
a program for further investigation. 



2. NUMERICAL TECHNIQUES 

2.1. Hydrodynamics 

For the 2D simulations presented in this paper, we use 
BETHE-hydro (Murphy & Burrows 2008a,b), an Arbitrary 
Lagrangian-Eulerian (ALE) hydrodynamics code. The ba- 
sic equations of hydrodynamics are the conservation of mass, 
momentum, and energy: 
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p is the mass density, v is the fluid velocity, $ is the gravi- 
tational potential, P is the isotropic pressure, e is the specific 
internal energy, and d/dt = d/dt + v • V is the Lagrangian time 
derivative. In this work, the neutrino heating, TL, and cooling, 
C, terms in eq. (3) are assumed to be 
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Note that these approximations for heating and cooling by 
neutrinos (Bethe & Wilson 1985; Janka 2001) depend upon 
local quantities and predefined parameters. They are p, tem- 
perature (T), the distance from the center (r), the electron- 
neutrino temperature (T Ve ), and the electron-neutrino lumi- 
nosity (L„ c ), which we give in units of 10 52 erg s" 1 . By us- 
ing eqs. (4) and (5), we gain considerable time savings by 
approximating the effects of detailed neutrino transport. For 
all simulations, we set T„ t = 4 MeV. In eq. (4), it has been as- 
sumed that L Ve = Ly e and that the mass fractions of protons and 
neutrons sum to one. Therefore, the sum of the electron- and 
anti-electron-neutrino luminosities is L u ,v, = 2L Ve . Closure for 
eqs. (1-3) is obtained with the finite-temperature nuclear EOS 
of Shen et al. (1998) to which we have added the effects of 
photons, electrons, and positrons. As such, the EOS has the 
following dependencies: 



P = P(p,e,Y e ), 



(6) 



where Y e is the electron fraction. Given this Y e dependence, 
we also solve the equation: 



dY^ 
dt 



(7) 



where r e is the net rate of Y e change. 

The heating and cooling functions, eqs. (4) & (5), are valid 
for optically thin regions only. To mimic the reduction in heat- 
ing and cooling at depth where the PNS is opaque to neutri- 
nos, we multiply by e~ Tdl . In ID simulations, the effective 

7 The neutrino distribution is entirely forward peaked and the flux factor is 
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optical depth for electron- and anti-electron-type neutrinos is 
calculated by 

/CO 
K e s(r)dr, (8) 

where the effective opacity, K e ff is given by eq. (15) or (16) of 
Janka (2001), which we reproduce without derivation: 

K*f) « 1.5 X 10-%, (^^3) (^)-- , 

(9) 

where X„ p , the composition weighting, is ~ j{Y n + Y p ) and T,, 
and Y p are the number fractions of free neutrons and protons. 
Because the radial density profile after bounce is monotonic 
and, to good approximation, a series of power-laws, r e ff is 
roughly a function of local density only. We use our ID sim- 
ulations to calibrate this r-p relationship and parametrize r as 
function of p for the 2D simulations. 

The heating and cooling terms are set to zero during col- 
lapse and are turned on after bounce. We have experimented 
with several methods to do this, but find that they give sim- 
ilar results provided that these terms "turn on" within a few 
milliseconds of bounce. For convenience, we include heating 
and cooling once the maximum value of r is greater than 100. 

Using BETHE-hydro (Murphy & Burrows 2008a), we 
solve eqs. (1-3) in two dimensions by the Arbitrary 
Lagrangian-Eulerian (ALE) method. To advance the discrete 
equations of hydrodynamics by one time step, ALE meth- 
ods generally use two operations, a Lagrangian hydrodynamic 
step followed by a remap. The structure of BETHE-hydro's 
hydrodynamic solver is designed for arbitrary-unstructured 
grids, and the remapping component offers control of the time 
evolution of the grid. Taken together, these features enable the 
use of time-dependent arbitrary grids to avoid some unwanted 
features of traditional grids. 

For the calculations presented in this paper, we use this flex- 
ibility to avoid the coordinate singularity of spherical grids in 
two dimensions. While spherical grids are generally useful 
for core-collapse simulations, the convergence of grid lines 
near the center places extreme constraints on the time step 
via the Courant-Friedrichs-Levy condition. A common rem- 
edy is to simulate the inner ^10 km in ID or to use an inner 
boundary condition. Another approach, which has been used 
in VULCAN/2D simulations (Livne 1993; Livne et al. 2004; 
Burrows et al. 2007c), is to avoid the singularity with a grid 
that is pseudo-Cartesian near the center and smoothly transi- 
tions to a spherical grid at larger radii. 

In this paper, we use a similar grid, the butterfly mesh (Mur- 
phy & Burrows 2008a,b) interior to 34 km and a spherical grid 
exterior to this radius. See Murphy & Burrows (2008a) for 
an example. For all 180° simulations, the innermost pseudo- 
rectangular region is 50 by 100 zones, and the region that 
transitions from Cartesian to spherical geometry has 50 ra- 
dial zones and 200 angular zones. The outermost spherical 
region has 200 angular zones and 500 radial zones which are 
logarithmically spaced between 34 km and 8000 km in most 
cases. Due to EOS constraints, the outer boundary goes out 
to 6000 km for the 12-M© model. The butterfly portion has 
an effective radial resolution of 0.34 km with the shortest cell 
edge being 0.28 km. In 2D simulations, we do not remap ev- 
ery time step. Since the time step is limited by the large sound 
speeds in the PNS core, we can afford to perform several La- 
grangian hydrodynamic solves before remapping back to the 



original grid. By remapping rarely, we gain considerable sav- 
ings in computational time. 

We make several approximations to expedite the calcula- 
tions. For one, gravity is calculated via g = —GM- mt /r 2 , where 
Mi n t is the mass interior to the radius r. Furthermore, since we 
do not solve the neutrino transport equations, we do not obtain 
electron/positron capture and emission rates, which are neces- 
sary for self-consistent Y e profiles and evolution. Instead, we 
use a Y e vs. p parametrization that is quantitatively accurate 
during collapse and gives the correct qualitative trends after 
bounce. Liebendtirfer (2005) observed that ID simulations 
including neutrino transport produce Y e values during collapse 
which are essentially a function of density alone. This allows 
for a parametrization of Y e as a function of p. To change Y e , we 
use results of ID SESAME (Burrows et al. 2000; Thompson 
et al. 2003) simulations to define the function Y e (p), and we 
employ the prescription of Liebendtirfer (2005) to calculate 
local values of r e . Though this prescription is most appro- 
priate during collapse, we also apply it after bounce to con- 
tinue deleptonizing postshock material at later times. While 
this is not entirely consistent with our heating and cooling pre- 
scription, the most important effects of deleptonization are in- 
cluded without the need for expensive, detailed neutrino trans- 
port. 

The approximations for neutrino-matter interactions, elec- 
tron capture rates, and gravity give simulations that reproduce 
the primary features of the core-collapse problem. These in- 
clude a PNS, a stalled shock at ^200 km, a gain radius at 
^100 km (that divides the gain region above with net heating 
from the cooling region below), prompt convection, PNS con- 
vection from ^20 to ^45 km, postshock convection, and the 
SASI. While the Y e -p parametrization is designed to give the 
proper Y e profile during collapse, we do not reproduce the Y e 
trough at later times that extends from ~40 km to the shock 
(e.g., Dessart et al. 2006). Dessart et al. (2006) conclude that 
the negative lepton gradient from ~15 to ^30 km in simula- 
tions employing consistent neutrino transport drives inner the 
PNS convection. In our simulations, PNS convection is driven 
by a negative entropy gradient. Although ID simulations that 
solve neutrino transport more consistently show a similar neg- 
ative entropy gradient, the depth and width are smaller in 
those simulations (compare Figs. 13 and 6 of Murphy & Bur- 
rows 2008b; Liebendtirfer et al. 2005). Furthermore, our PNS 
convection velocities are roughly a factor of two larger than 
those obtained by Dessart et al. (2006). In summary, while 
we obtain a PNS convective region, our PNS convection is 
predominantly driven by an entropy gradient rather than by a 
negative lepton gradient and is more extended and vigorous 
compared to more realistic calculations. However, we show 
in later sections that though the GW signal due to PNS con- 
vection is important, it is not an overwhelmingly dominant 
feature of GW emission from CCSNe. Therefore, the differ- 
ences in PNS convection should not change the qualitative 
conclusions of this paper. 

2.2. Progenitor Models 

As initial conditions, we use the 12-, 15-, 20-, and 40- 
M Q core-collapse progenitor models of Woosley & Heger 
(2007), where these masses correspond to the zero-age-main- 
sequence (ZAMS) masses. In Fig. 1, we plot the mass ac- 
cretion rate vs. time after bounce (solid black lines) for each 
of these models. In general, because the more massive pro- 
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genitors have extended core structures, the accretion rates for 
them are higher at a given time. Since higher accretion rates 
require larger neutrino luminosities for successful explosions 
(Murphy & Burrows 2008b), this translates into higher crit- 
ical neutrino luminosities for the more massive progenitors. 
In fact, if we assume that only the neutrino mechanism leads 
to explosion, the required neutrino luminosity for successful 
explosions of the 40-M Q progenitor is so high, > 10 53 erg s" 1 , 
we doubt that it will explode at all and suspect that the non- 
explosive simulation using L Ve = 6 x 10 52 erg s" 1 represents 
the likely outcome of core collapse for the 40-M Q progenitor. 
In general, these trends with progenitor mass will translate 
to trends in the GW signal which we discuss in forthcoming 
sections. 

Although core structure follows a general trend with the 
progenitor's ZAMS mass, the correlation is known to be 
somewhat chaotic (Woosley et al. 2002). For example, a sub- 
stantial drop in the 20-M Q model's accretion rate correspond- 
ing to the interface between the Si and O burning shells tem- 
porarily reduces the accretion below even the 12- and 15-M Q 
progenitors at ^300 ms past bounce. This has consequences 
for the development of the SASI, explosion, and the GW sig- 
nal. 

3. RESULTS 

3.1. Gravitational Wave Extraction and Analysis 

We extract the gravitational wave signal from our hydrody- 
namic simulations via the slow-motion, weak-field formalism 
(e.g., Misner et al. 1973) and consider the dominant mass- 
quadrupole only. In this approximation, the gravitational 
wave field is directly proportional to the transverse-traceless 
(TT) part of the second time derivative of the reduced mass- 
quadrupole tensor^. Specifically, 
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In the above, D is the observer distance to the source of emis- 
sion. All other variables and constants have their usual mean- 
ing and the Einstein sum convention is assumed. For numer- 
ical convenience, we employ the so-called first-moment-of- 
momentum-divergence (FMD) recast of eq. (10) proposed by 
Finn & Evans (1990), which employs the continuity equation 
to remove one time derivative. In axisymmetry, 4-jj reduces 
to one independent component (4- zz ) and, adapted to spheri- 
cal coordinates, the FMD formula for its first time derivative 
reads 
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where P2(co&8) is the second Legendre polynomial in cosf? 
and v r and vg are the fluid velocities in the radial and lateral 
direction, respectively. The axisymmetric 8 gravitational wave 



strain h+ = hgg (dropping the TT superscript) is then given by 
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where a is the angle between the symmetry axis and the line 
of sight of the observer and the second time derivative is ob- 
tained via straightforward differentiation 9 . For our simula- 
tions, this approach yields good results and we do not find it 
necessary to employ the modification of eq. (10) proposed by 
Blanchet et al. (1990) that removes both time derivatives, but 
introduces a sensitive dependence on derivatives of the New- 
tonian gravitational potential, which we treat only in spherical 
fashion. 

In our analysis of the GW signature, we not only consider 
the dimensionless GW strain /?+, but also compute the total 
emitted energy in GWs, given by 



F - 3 ° 

£gw -io? 



d' 



dt . 



and the GW spectral energy density via 



dE(yw 3 G 



df 



--{2KfY\A\ 



(14) 



(15) 



where A denotes the Fourier transform of A = 4t^- computed 



via 



A(f)= / A(t)e- 2 ^'dt 



(16) 



At this point, it is useful to define for future reference the di- 
mensionless characteristic gravitational wave strain (Flanagan 
& Hughes 1998), in terms of the GW spectral energy density, 



2 G 1 c/£gw 



*char " 



r 2 „3 



c 3 D 2 df 



(17) 



For signals with relatively stable frequencies and ampli- 
tudes, Fourier transforms and their energy spectra are ade- 
quate frequency analysis tools. However, for signals with 
time-varying amplitudes and frequencies, a short-time Fourier 
transform (STFT) is more appropriate. The STFT of A(t) is 



S(f,T): 



A(t)H(t-T)e- 2nif 'dt , 



where r is the time offset of the window function, H(t 
We use the Hann window function: 



W-T): 



1 +COS 



' n(t-r) 
, St 



)) 



for 
for 



■t|<! 

•r|>| 



(18) 
-r). 

(19) 



where St is the width of the window function. The analog of 
the energy spectrum of the Fourier transform is the spectro- 
gram, |5(/,r)| 2 . Using the spectrogram, we define an analog 
to the energy emission per frequency interval (eq. 15): 



d -^L(f,r)= 3 -^(2n f y-\S(f,r)\ 2 . 



(20) 



We emphasize that the GW strains reported in this paper are 
based upon matter motions alone and do not include the low 



8 Note that in axisymmetry h x = everywhere and h+ is non-zero only 
away from the axis of symmetry of the system. 



9 Note that due to a misprint in Eq. (38) of Finn & Evans 1990 their ex- 
pression is a factor of 2 smaller than our expression in Eq. 13. 
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frequency signal that results from asymmetric neutrino emis- 
sion (Burrows & Hayes 1996; Miiller & Janka 1997). Ac- 
curate calculations of asymmetric neutrino emission require 
multi-D, multi-angle neutrino transport to capture the true 
asymmetry of the neutrino radiation field (see, e.g., Ott et al. 

2008) . Our choice to parametrize the effects of neutrino trans- 
port by local heating and cooling algorithms is based upon as- 
sumptions of transparency, which ignore diffusive effects and 
would exaggerate the asymmetries and resulting GWs. For 
example, Kotake et al. (2007) estimated the neutrino GW sig- 
nal using a similar heating and cooling parametrization and 
obtained GW strain amplitudes that are ~100 times the mat- 
ter GW signal. However, with an improved ray-tracing-based 
method, the same authors find much smaller amplitudes that 
are larger than those due to matter motions by only a factor 
of a few (Kotake et al. 2009). This is in agreement with the 
GW estimates of Marek et al. (2009) who used ID ray -by-ray 
neutrino transport and coupled neighboring rays in 2D hydro- 
dynamic simulations. 

Studying the matter GW signal alone is worthwhile. Al- 
though the neutrino GW strain amplitudes can be as large or 
even larger than the contribution by matter (Burrows & Hayes 
1996; Miiller & Janka 1997; Miiller et al. 2004; Marek et al. 

2009) , the typical frequencies, /, of the neutrino GW sig- 
nal (~10 Hz or less) are typically much lower than the fre- 
quencies of the matter signal (>100 Hz). Consequently, the 
GW power emitted, which is proportional to / 2 , can be much 
higher for the matter GW signal. Furthermore, although fu- 
ture GW detectors (e.g. Advanced LIGO) will have improved 
sensitivity at low frequencies, current detectors have response 
curves that are not sensitive to the lower frequencies of the 
neutrino GW signal. 

3.2. Signatures in the GW Strain 

In Fig. 1, we plot the GW strain (eq. 13) times the distance 
to a 10 kpc source, h+D, vs. time after bounce for all sim- 
ulations. Though there is some diversity in amplitude and 
timescale among these GW strains, there are several recur- 
ring features that exhibit systematic trends with mass and neu- 
trino luminosity. We illustrate these features in Fig. 2 with 
the GW strain of the simulation using the 15-M Q progenitor 
and L Ve = 3.7 x 10 52 erg s . Before bounce, spherical col- 
lapse results in zero GW strain. Just after bounce the prompt 
shock loses energy and stalls, leaving a negative entropy gra- 
dient that is unstable to convection. Because the speeds of 
this prompt convection are larger than those of steady-state 
postshock or PNS convection afterward, the GW strain ampli- 
tude rises to h+D ~ 5 cm during prompt convection and settles 
down to ~1 cm roughly 50 ms later, which is consistent with 
the results of Ott (2009b) and Marek et al. (2009). Later in 
this section, we show that during both phases, convective mo- 
tions in postshock convection above the neutrinosphere and 
PNS convection below it contribute to the GW strain. Since 
nonlinear SASI oscillation amplitudes increase around 550 
ms past bounce, the GW signal strengthens from h+D ~ 1 to 
10 cm and is punctuated by spikes that are coincident in time 
with narrow plumes striking the PNS "surface" (at ^50 km). 
Marek et al. (2009) also noted this correlation. 

The final feature after ^800 ms is associated with explo- 
sion. The signatures of explosion are twofold. First, dur- 
ing explosion, postshock convection and the SASI subside 
in strength and the higher frequency (^300-400 Hz) oscil- 



lations in h+D diminish. Second, global asymmetries in mass 
ejection result in long-term and large deviations of the GW 
strain. In Fig. 2, a monotonic rise of h+D to nonzero, specif- 
ically positive, values corresponds to a prolate explosion in 
this simulation. This is similar to the "memory" in the GW 
signal of asymmetric neutrino emission (Burrows & Hayes 
1996; Miiller & Janka 1997). When the explosion is spher- 
ical, the strain drops to zero and remains there, and when it 
is oblate, the strain maintains negative values. Examples of 
h+D curves showing prolate, oblate, and spherical explosions 
are shown in Fig. 3. The simulation using a 20-M Q progenitor 
and L Ve = 3 .4 x 1 52 erg s" 1 (gray line) exploded with an oblate 
structure, the simulation with M = 12 M Q and L Ve = 3.2 x 10 52 
erg s -1 (light gray line) has a prolate explosion, and the simu- 
lation with M = 12 M Q and L Ve = 2.2 x 10 52 erg s" 1 explodes 
almost spherically. 

Figure 4 shows snapshots of the entropy (in units of 
fc^/baryon) distribution during explosion for the three mod- 
els highlighted in Fig. 3. Lighter shades (warmer colors in 
the online version) represent higher entropies, while darker 
shades (cooler colors in the online version) represent lower 
entropies. The general shapes of the matter interior to the 
shocks is oblate (M = 20 M Q and L Ue = 3.4 x 10 52 erg s" 1 ), 
prolate (M = 12 M Q and L Vt = 3.2 x 10 52 erg s" 1 ), and spheri- 
cal (M = 12 M Q , L Ve = 2.2 x 10 52 erg s" 1 ; and thin black line). 
We emphasize that the GW signal is sensitive to £ = 2 acceler- 
ations of matter and is somewhat blind to differences in com- 
position or higher order asymmetries. Therefore, even though 
the GW signal may indicate that the explosion is in general 
"spherical," "oblate," or "prolate," the entropy (hence temper- 
ature and composition) distributions may not be. 

In Figs. 5 and 6, we localize the source of GW emission. 
Figure 5 shows the spatial distribution of 



1 d 

P 2 (cos B) v r + - -t^Pi (cos 9) v e 



(21) 



which is the time derivative of the integrand of eq. (12) and 
determines the GW strain (eq. 13). As in Fig. 2, these re- 
sults are for the 15-M Q progenitor and L Ve = 3.7 x 10 52 erg 
s" 1 . Lighter shades (brighter colors in the online version) cor- 
respond to a larger contribution to the GW strain integral. The 
left panel is at 615 ms past bounce, and represents the phase 
with a stalled shock (~250 km) and vigorous postshock con- 
vection/SASI motions. This is our first indication that motions 
below the gain radius (^100 km), in particular deceleration of 
plumes and PNS convection, are the strongest sources of GWs 
before explosion. 

The right panel is at 1060 ms past bounce and shows that ac- 
celerations at an aspherical shock (~4000 km at this time) are 
responsible for the "memory" signature at late times. Using 
the appropriate shock velocities and asymmetries and post- 
shock densities and velocities in eqs. (12) and (13), we obtain 
the correct order-of-magnitude for the memory signal. How- 
ever, during explosion a high entropy wind emerges that en- 
counters the swept-up material of the primary shock and pro- 
duces secondary shocks. We find that these secondary shocks 
and their asymmetric structure can add a non-trivial contribu- 
tion to the "memory" signal. The strength of these winds is a 
strong function of the neutrino luminosity (Qian & Woosley 
1996; Thompson et al. 2001). Since we employ a constant 
neutrino luminosity, an accurate characterization of the late- 
time "memory" (such as saturation) are beyond the scope of 
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this paper. 

Figure 6 quantifies the sources of GW radiation and their 
spatial distribution as a function of radius. The model shown 
is the same as in Figs. 2 and 5. For reference, the entire GW 
signal is shown in both panels (solid-black line). The signals 
originating from > 50 km (orange, online version) and < 50 
km (blue, online version) are shown in the top panel. Below 
50 km, PNS convection dominates the motions, and above 
this radius postshock convection and the SASI dominate. As 
expected, most of the signal associated with prompt convec- 
tion originates in the outer convective zone, though motions 
below 50 km and presumably from PNS convection account 
for a fair fraction (^20%). Afterward, the contributions to 
the GW amplitude from below and above 50 km are compara- 
ble. This suggests that motions associated with both PNS con- 
vection and the postshock-convection/S ASI region contribute 
significantly to the GW signal. The "memory" signature of 
the GW strain during explosion clearly originates from the 
outer (exploding) regions. Interestingly, once the model ex- 
plodes, and the nonlinear postshock convection/SASI motions 
subside, the GW signal from below 50 km diminishes as well, 
though PNS convective motions do not. 

To further understand the origin of GW emission, the bot- 
tom panel of Fig. 6 refines the spatial origin of the GW emis- 
sion into five regions. We plot the GW strain from five over- 
lapping regions, and each extends from the center to different 
outer radii (30, 40, 50, 60, and 100 km). PNS convection ex- 
tends from ^20 to ^40 km and the gain radius is ^100 km. 
However, as we explain in section §3.4, the turbulent motions 
of postshock convection/SASI penetrate below the gain ra- 
dius to radii of ~ 60 km during the most vigorous phases. 
Therefore, the partial GW signals with outer radii of 30 and 
40 km contain signals only due to PNS convection; 50 km 
encompasses PNS convection and gravity waves that are ex- 
cited by the overlying convection/SASI; 60 km encompasses 
PNS convection, gravity waves, and the deepest penetration of 
the convection/SASI plumes (see §3.4); and 100 km encom- 
passes all of these contributions and the gain radius. First, 
we note that extending the outer radius from 60 to 100 km, 
the gain radius, adds very little signal during the most vig- 
orous postshock convection/SASI phase from 550 to 800 ms 
past bounce. Therefore, the most relevant motions are those 
at ~60 km and below, which include the postshock convec- 
tion/SASI plume decelerations, excited gravity waves, and 
PNS convection. Furthermore, the strength of the GW sig- 
nal from all radii increases and decreases in synchrony with 
the vigor of postshock convection/SASI motions, which are 
restricted to radii above ~50 km. This suggests that the influ- 
ence of the convection/SASI diminishes only gradually with 
depth. Even though the mechanism for convection/SASI and 
their nonlinear motions are above 50 km, they influence PNS 
convective motions below 50 km, which results in GW emis- 
sion whose vigor is ultimately due to the postshock convec- 
tion/SASI motions. 

Hence, we conclude that features of the SASI and post- 
shock convection, in particular the plumes striking the PNS 
"surface", are ultimately responsible for the strong GW sig- 
nal from ^550 ms until explosion at ^800 ms after bounce. 
Some of the signal comes from the accelerations at this inter- 
face. The rest comes from gravity waves and motions within 
the PNS (< 50 km) that are excited by the plumes striking this 
layer. Whether these motions are excited g-waves, enhanced 



PNS convection, or a combination of both is difficult to dis- 
tinguish in our simulations. A focused investigation on this 
aspect is warranted. 

Miiller et al. (2004) conclude that the GW signal from PNS 
convection is a few factors smaller in amplitude than the GWs 
from postshock convection. At first sight, this conclusion ap- 
pears at odds with our results that the signals above and below 
50 km are roughly equal. However, to calculate the GW sig- 
nal from PNS convection, Mtiller et al. (2004) used results 
from a simulation that extended out to only 60 km (Keil et al. 
1996), which ignores any influence that postshock-convective 
or SASI motions have on PNS convection. This is consis- 
tent with the GW signal we obtain after explosion, when post- 
shock convection and SASI go away, but PNS convection re- 
mains. These results emphasize the importance of calculating 
PNS convection in the full context of core-collapse dynamics. 

We now return to Fig. 1 and analyze the trends with progen- 
itor mass and L„,. Each panel represents a single progenitor 
model and the curves are labeled by L Vc . The accretion rates 
(in Mq s" 1 ) at the stalled shock of non-exploding models are 
shown for comparison. The prominent features described in 
Fig. 2 are apparent in all models, except the run using M = 40 
M and L Ue = 6.0 x 10 52 erg s" 1 which does not explode and 
manifests the larger amplitudes associated with the nonlinear 
SASI at much later times. Otherwise, lower-luminosity sim- 
ulations take longer to develop the nonlinear-SASI GW sig- 
nal and explode at later times, but all runs seem to saturate 
at roughly similar amplitudes. This is consistent with the re- 
sults of Murphy & Burrows (2008b), in which they note that 
lower luminosity runs take longer to reach saturated shock 
oscillations, but that the saturated amplitudes are similar for 
the different neutrino luminosities. Similarly, Kotake et al. 
(2009) conclude that the GW amplitudes are independent of 
neutrino luminosity. Finally, we note that the significant drop 
in accretion rate in simulations using the 2O-M progenitor, 
instigates strong nonlinear SASI motions and associated GW 
signatures. This might suggest using GWs as a diagnostic 
for the location of shell burning or compositional interfaces. 
However, very few of the other simulations show such a con- 
vincing correlation. 

The neutrino luminosities required to explode the 4O-M 
progenitor model are much larger than simulations obtain that 
consistently calculate neutrino transport (Dessart et al. 2006; 
Ott et al. 2008; Marek & Janka 2009). Therefore, if other 
mechanisms fail to explode higher mass stars, we suspect that 
the simulation using L Ve = 6 x 10 52 erg s _1 , which does not 
explode by the neutrino mechanism, gives an upper limit on 
the most likely GW emission strength for massive progenitors 
that collapse to form black holes without explosion. Taken at 
face value, the low convective and SASI motions of this model 
imply that, in the absence of rapid rotation, the GW power 
emitted by the most massive progenitors that form black holes 
could be quite low. 

3.3. GW Energy Spectra and Detection 

If rotation rates are low and PNS core g-modes are weak, 
then postshock-convection/SASI motions are the most prob- 
able source of GWs. In this scenario, the GW emission of 
Galactic CCSNe (^10 kpc) is probably too weak to detect 
the GW strain waveform directly as a time series. Instead, 
GW observers will search for excess spectral power above 
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the detector noise in time-frequency maps (e.g., Abbott et 
al. 2009). The energy spectra vs. frequency for all simula- 
tions is plotted in Fig. 7 assuming a distance to the source, D, 
of 10 kpc. Rather than showing dEcw/df, we plot /i c h ar (eq. 
17), which is proportional to y/dEcw/df. As in Fig. 1, each 
panel shows the spectra from simulations for the same progen- 
itor model, and within each panel the spectra are color-coded 
and labeled by L Vr . For reference, the noise thresholds for ini- 
tial LIGO (solid black line, Gustafson et al. 1999), Enhanced 
LIGO (dot-dashed-black line, Adhikari 2009), and Advanced 
LIGO (burst mode, dashed-black line, Shoemaker 2006) are 
included. 

Characteristically, the spectra show broad peaks with max- 
imum power at ~300 to ~400 Hz. The lower and higher fre- 
quencies correspond to lower and higher progenitors masses, 
respectively. Later in this section and §3.4, in analyzing the 
spectra as a function of post-bounce time, we show in Fig. 
8 that simulations with lower L Vt explode later and obtain 
higher frequencies at later times. However, the spectra in Fig. 
7, which are calculated using the entire time sequence, do not 
clearly show such a correlation. The spectra associated with 
12-, 15-, and 20-M Q progenitors show a second, weaker peak 
at ^100 Hz that we associate with prompt convection. Gen- 
erally, all spectra are only marginally above or just below the 
design noise threshold for Initial LIGO. On the other hand, 
all simulations have power above the design noise threshold 
for Enhanced and Advanced LIGOs, with maximum spectral 
signal-to-noise ratios of ~5 and ~ 10-20, respectively. 

Characteristic GW amplitudes, frequencies, optimal theo- 
retical single-detector signal-to-noise ratios, and energies for 
all simulations presented in this paper are listed in Table 1. 
The first two columns identify the progenitor mass and neu- 
trino luminosity of the simulation. Following are the maxi- 
mum amplitude of the GW strain (|/j+, rms |), the signal-to-noise 
ratios with respect to the Initial, Enhanced, and Advanced 
LIGO sensitivity curves (S/N L tgo, S/N eL iGo, and S/N a d v uGo)- 
We also provide values for the emitted GW energy (£bw) and 
h c hm (eq- 17) for the maxima of the characteristic strain spec- 
tra and the frequencies / c har,max at which they are located. The 
values in this table do not present clean, monotonic trends 
with progenitor mass or neutrino luminosity. However, we 
note some general trends. More massive progenitors tend to 
produce higher frequencies. 10 On the other hand, for a given 
progenitor model, higher neutrino luminosities give lower fre- 
quencies. In §3.4, we explain that the characteristic frequency 
of GWs is connected with the buoyancy frequency in the post- 
shock region, which is higher for more massive progenitors, 
increases with time, and is largely independent of neutrino lu- 
minosity. This explains the first correlation with progenitor 
mass, but to explain the anti-correlation with neutrino lumi- 
nosity, we note that the higher luminosity simulations explode 
earlier when the buoyancy frequency is lower. 

In general, the typical GW strain for a source at 10 kpc is 
a few times 10~ 22 to 10~ 21 , and typical values of £gw range 
from a few times 10~ n to a few times 10~ 10 M Q c 2 . Ko- 
take et al. (2009) conclude that the stochastic properties of 
the SASI preclude any relationship between the neutrino lu- 
minosity and the GW amplitude. While there appears to be 
no trends, we note several competing dependencies on neu- 

10 Burrows et al. (2007c) show in Table 1 of their paper a similar trend in 
that higher progenitor masses produce higher shock-oscillation frequencies, 
though at much smaller values. 



trino luminosity that result in a non-monotonic and compli- 
cated relationship. In §3.4, we connect the GW amplitude to 
the speed of plumes, v p , at the base of the SASI region. In 
addition, we note that although lower luminosity simulations 
take longer to saturate the SASI motions, the saturation value 
eventually achieved is similar for all runs. Since most of the 
GW power is emitted during the phase of strongest postshock 
convection/SASI motions, the total energy emitted is deter- 
mined by the duration of this phase, which starts earlier for the 
higher neutrino luminosity, but also ends much earlier at the 
higher luminosities, so that the total duration of the stronger 
SASI phase, and therefore, £gw is a non-monotonic function 
of neutrino luminosity. 

Focusing on the energy spectrum for an entire GW strain 
data set is most appropriate when the signal is regular in 
both amplitude and frequency. However, Figs. 1 and 2 show 
that the amplitude and frequency change substantially in these 
simulations. In this case, a time-frequency analysis, or spec- 
trogram, is more appropriate (see §3.1). To obtain energy 
spectra vs. time we take the Fourier transform of the GW 
strain convolved with a Hann window function with width 
^20 ms and sample the time domain, r in eq. (18), at intervals 
of 2 ms. In Fig. 8, we show color maps of dE^ w /df (eq. 20) 
vs. frequency and time. The warmer colors in Fig. 8 reflect 
higher power, and vice versa for cooler colors. In general, 
simulations that explode later (i.e. with higher mass progen- 
itors and/or lower neutrino luminosities) have GW power at 
higher frequencies. During explosion, the frequency at peak 
power drops to 10s of Hz, owing to asymmetric mass ejection. 

Unlike in the spectra, which are calculated using the en- 
tire time domain (Fig. 7), prompt convection, nonlinear SASI, 
SASI plumes, and explosion are quite apparent in these time- 
frequency plots. Figure 9 labels these features in the spectro- 
gram for the simulation using M = 15 M© and L Ve = 3.7 x 10 52 
erg s -1 . The three features showing the largest power are 
prompt convection just after bounce, nonlinear SASI mo- 
tions, and explosion. Though the power declines between 
prompt convection and the start of the nonlinear SASI, the 
frequency at peak power increases monotonically from ^100 
Hz at bounce to ~300-400 Hz just before explosion. In the 
next section, we propose that the characteristic timescale for 
the deceleration of plumes striking the PNS "surface" deter- 
mines the GW peak frequencies. The solid black line in Fig. 9 
shows our analytic estimate for this characteristic frequency, 
f p , and that it agrees with the frequency and time evolution of 
the GW signal from simulation. 

3.4. Source of GWs and a Theory for Their Characteristic 
Frequencies and Amplitudes 

Figure 6 shows that motions associated with the broad re- 
gion which includes PNS convection, postshock convection, 
and the SASI contribute to the overall GW signal. Marek et al. 
(2009) come to a similar conclusion. However, in § 3.2 and 
§ 3.3, we noted that the strength of the GW signal waxes and 
wanes as the vigor of postshock convection and SASI mo- 
tions rises and falls. Moreover, spikes in the GW strain often 
correlate in time with downdrafts (plumes) striking the PNS 
"surface" (also noted by Marek et al. 2009). Hence, we ar- 
gue that the strength and characteristic frequencies of the GW 
emission are determined predominantly by the deceleration of 
the plumes striking the PNS "surface." 

In other contexts, the "surface" of the PNS is defined as the 
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position of the neutrinosphere, the location of a steepening 
density gradient, or the narrow region over which the aver- 
age radial velocity of the accreting matter approaches zero. 
These definitions roughly coincide in space and are in fact 
intimately related, but they neither explain nor quantify the 
strong decelerations of the downdrafts in an otherwise contin- 
uous medium. We note, however, that buoyancy forces and 
their radial profile at the boundary between convective and 
stably-stratified regions explain the deceleration of the plumes 
and, hence, the GW amplitudes and frequencies. We, there- 
fore, define the PNS "surface" as the region where the nega- 
tively buoyant downdrafts become positively buoyant and re- 
verse their downward motions. 

An important quantity in analyzing convection and buoy- 
ancy is the square of the Brunt-Vaisala (or buoyancy) fre- 
quency 

2 ( GM r \ ( 1 dlnP dlnp 



T] d\nr d\nr 



(22) 



where M r is the mass interior to the radius, r, and the thermo- 
dynamic derivative, Ti = (dlnP/dlnp)s, is evaluated at con- 
stant entropy, S. In the local dispersion relation for waves, this 
corresponds to the characteristic frequency associated with 
gravity waves (not to be confused with gravitational waves), 
whose dominant restoring force is buoyancy, which sets the 
frequency scale. Alternatively, N 2 <Q indicates linear insta- 
bility and is in fact the Ledoux (and Rayleigh-Taylor) con- 
dition for convection. Naively, one might expect the radius 
where A^ 2 changes sign from negative to positive to mark the 
lower boundary of the convection zone. However, the Ledoux 
condition is derived from a linear analysis and ignores im- 
portant nonlinear effects. In practice, convective plumes have 
momentum when they reach this boundary and penetrate be- 
yond the Ledoux condition boundaries. This is called over- 
shoot and has a long history and many proposed prescriptions, 
though many of these are either inappropriate for the dynam- 
ics in core-collapse SNe or lack adequate physical motivation, 
resulting in a free parameter (see Meakin & Arnett (2007b), 
Arnett et al. (2009), and references therein). For this paper, we 
adopt an analysis based upon buoyancy and the bulk Richard- 
son number (Meakin & Arnett 2007b), which is physically 
motivated and contains no free parameters. 

Beyond the boundaries defined by the Ledoux condition, A^ 2 
is positive and, therefore, coherent plumes that penetrate this 
boundary experience a buoyancy force back toward the cen- 
tral regions of convection. In this context, the buoyancy ac- 
celeration, b{r), felt by a Lagrangian mass element displaced 
from r,- to r is related to the buoyancy frequency by 



b(r)= / N 2 (r)dr. 



(23) 



For plumes approaching the lower boundary of postshock 
convection, the magnitude of the integrand is largest where 
Af 2 > 0. Therefore, we take as r, the lower Ledoux condition 
boundary. With a plume velocity in the radial direction, v p , 
at r, as initial conditions and eq. (23) as the position depen- 
dent acceleration, we integrate the plume's equation of mo- 
tion to determine the depth of penetration, D p . Because the 
sources of gravitational waves are time changing quadrupole 
motions, the inverse of the characteristic time of the buoy- 
ancy impulse for each plume should correspond roughly to the 
peak frequencies in the GW spectra and spectrograms (Figs. 
7 and 8). This timescale, t p , is defined as the HWHM of the 



buoyancy acceleration pulse and the associated frequency is 
fp = 1 /{2ntp) n . It is this definition of for which we give 
quantitative results in this paper. 

Using the bulk Richardson number, a dimensionless mea- 
sure of the boundary layer stiffness, we obtain an analytic de- 
scription of / ; , that is independent of, and corroborates, our 
method for calculating f p . There is a long tradition of us- 
ing this dimensionless parameter in atmospheric sciences to 
accurately describe the boundary layer in atmospheric con- 
vection (Fernando 1991; Fernando & Hunt 1997; Bretherton 
et al. 1999; Stevens 2002). More recently Meakin & Arnett 
(2007b) have successfully used similar arguments to describe 
the boundaries of convection in stellar evolution. The bulk 
Richardson number is 



Rh = 



AbD p 

V 



(24) 



where Ab is the change in buoyancy acceleration across the 
boundary, D p is a length scale such as the penetration depth, 
and v 2 is the typical velocity of the plume. This is not to 
be confused with the gradient Richardson number, R g , which 
is derived from a linear analysis and is a condition for shear 
instabilities in a stratified medium. While R g = 1/4 is the de- 
rived value that separates stability from instability in a linear 
analysis, Rb is a rough measure of the stiffness of a stable layer 
next to a convective layer in nonlinear flows. Approximately, 
it is the ratio of work per mass done by buoyancy (AbD p ) 
and the kinetic energy per mass of the plume (v 2 ). Therefore, 
Rb < 1 offers little resistance to penetration, Rb> 1 represents 
a stiff boundary, and where Rb ~ 1, the work done by buoy- 
ancy balances the plume's kinetic energy, causing it to turn 
around. Hence, Rb is appropriate in characterizing the pene- 
tration of downdrafts into the underlying stable layers. 

To estimate scales and trends, we approximate the integral 
in eq. (23), giving 



Ab ~ N l D 



p • 



(25) 



and we assume that R/, ~ 1, which we argued above is where 
buoyancy balances the plume's kinetic energy and gives the 
turn-around point. This gives (ND p /v p ) 2 ~ 1, or 



D„ 



N 



(26) 



We note that this scaling is similar to the gravity wave dis- 
placement amplitude derived by Meakin & Arnett (2007a), 
which they derive by equating the gravity wave pressure fluc- 
tuations to the ram pressure of the convective eddies at the 
boundary. An estimate of the characteristic timescale is given 
by the velocity divided by the acceleration, t p ~ v p / Ab. Sub- 
stituting in eqs. (25) and (26) gives 12 



Z7T 



(27) 



Interestingly, this characteristic frequency is insensitive to the 
plume velocity and penetration depth. Rather, f p is most sen- 
sitive to the buoyancy frequency at the turning point (AWn)- 

Although the above analysis suggests that the character- 
istic frequency is independent of \> p , the amplitude of the 

1 1 The factor of 2n arises because the definition of t p is similar to the 
a, in exp[— ? 2 /(2of)], and the Fourier transform of this Gaussian is another 
Gaussian with width at = l/(2nai). 

12 see footnote 11 
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GW strain is directly dependent upon these velocities. Very 
roughly, the GW strain is 

h + ~^pf p v p r 3 ArA[i, (28) 

where we use eqs. (12) and (13), Ap is the finite difference 
approximation of d cos 9, and we assume that d(pv p ) /dt ~ 
pdvp/dt and that dv p /dt ~ f p v p . Figure 10 compares h+D 
with the maximum downward plume speed below 120 km, 
v p , as a function of time after bounce. The spikes in v p cor- 
respond to the strongest plumes that strike the PNS "surface," 
and the baseline, below which v p does not dip, shows the time- 
averaged accretion of material onto the PNS. As the PNS ac- 
cumulates more mass, exerting stronger gravitational forces, 
the background and plume velocities evolve upward slightly 
until ^550 ms after bounce. At this time, the motions in 
the postshock convection/SASI region increase dramatically, 
leading to a significant rise in v p . After the start of explosion 
(>800 ms), it no longer makes sense to define a plume ve- 
locity, since the GW signal during explosion is dominated by 
explosion dynamics. Hence, we do not show v p after ~800 
ms. 

Initially, prompt convection dominates the signal and v p is 
not relevant. As prompt convection settles to steady-state con- 
vection, the plumes become more relevant. From ~150 ms to 
^550 ms, both v p and the GW amplitude are low, but as v p 
grows so does the GW amplitude. Strikingly, the increase in 
GW wave strain at ^600 ms coincides with the rise in v p . 
In fact, many, but not all, spikes in h+D coincide with rapid 
changes in v p . Furthermore after ^600 ms, substituting val- 
ues for typical flows near the boundary (e.g. p ~ 10" g cm -3 , 
v p ~ 4 x 10 9 cm s _1 , r ~ 60 km, Ar~D p ~ 20 km, f p ~ 300 
Hz, and Ap ~ D p /(irr)), into eq. (28) gives h+D ~ 5 cm, 
which is the order of magnitude of the GW strain. 

Figures 1 1 and 12 depict the buoyancy frequency, N (solid 
blue line, online version), (GM r / P) 1 ! 2 (dashed purple line, 
online version), and the radial velocity, v r (black points), as 
a function of radius at 300 and 620 ms after bounce for the 
2D simulation that uses the 15-M Q progenitor model and 
L Ve = 3.7 x 10 52 erg s" 1 . The velocities clearly show (from 
large to small radii), spherical infall, the asymmetric shock, 
the postshock-convection/SASI region, a stable layer, and 
PNS convection. At 300 ms (Fig. 1 1), the dynamics are dom- 
inated by postshock convection and mild SASI oscillations, 
while at 620 ms, the SASI oscillations have increased sub- 
stantially. Roughly, the buoyancy frequency scales with the 
average density, eq. (22). However, at the core, the spatial 
derivatives, and therefore N, approach zero. In the regions 
of PNS convection (~20 to ^40 km) and postshock convec- 
tion (>100 km), N 2 is negative. Green line segments in both 
figures indicate the regions that are unstable to convection by 
the Ledoux condition in ID simulations, and the penetration 
depths, D p , which are computed by integrating the plume's 
equation of motion, are labeled and indicated by an orange 
line. At 300 and 620 ms after bounce, D p at the base of the 
region of postshock convection is ~25 and ^39 km, respec- 
tively. These analytic penetration depths agree with those ob- 
tained from the numerical simulations. Note that the bound- 
aries of PNS convection in the 2D simulations coincide with 
the boundaries given by the Ledoux condition. This is due 
to the fact that the relatively large buoyancy frequency and 
small convective speeds imply very small overshoot depths 
(eq. 26) at the edges of PNS convection. The plume frequen- 



cies, fp, obtained from the HWHM timescale of the accelera- 
tion pulses, are ~196 Hz at 300 ms and ~272 Hz at 620 ms. 

To better understand the contributions to f p , we plot in 
Fig. 13 the characteristic frequency (/ ; ,), the buoyancy fre- 
quency at Dp (AWn/(27r)), and (GM r / r 3 ) 1/2 /(2tt) at D p as a 
function of time after bounce. Strikingly, and N tm - n / (2tt) 
agree quantitatively. Since is calculated using the equation 
of motion and is independent of the approximations and as- 
sumptions that led to eq. (27), the fact that is nearly equal 
to N tum /(2ir) for all times provides implicit confirmation of 
our approximations and assumptions in deriving the analytic 
result. Furthermore, we calculated the ratio N tum D p /v p (the 
square root of R/, after inserting the approximation that Ab ~ 
N 2 D p ) for all models (Fig. 14), found that it's value is around 
1 .5 to 2, and that it changes less than 30% during the simula- 
tion. This validates both the assumption that R/, ~ 1 charac- 
terizes where the plumes turn around and the approximation 
that Ab ~ N 2 D P . 

Initially, all three measures of frequency in Fig. 13 coin- 
cide. However, while N ium /{2ir) tracks the time-evolution of 
f p , (GM r /r 3 y/ 2 /(2n) diverges from the other two measures. 
By plotting both N lvm and (GM r /r 3 )^ 2 , which is the inverse 
free-fall time, we highlight global (inverse free-fall time) and 
local (difference of logarithmic derivatives) properties of the 
PNS that contribute to eq. (22). The inverse free-fall time is 
proportional to the square root of the average density (com- 
pactness), which is strongly dependent upon the dense-matter 
EOS. Marek et al. (2009) reported a change in peak GW fre- 
quency with stiff and soft dense-matter equations of state, 
though they make reference to only compact PNSs without 
explaining the timescales in the context of the buoyancy fre- 
quency. Fig. 13 shows that both the global compactness and 
the local gradients at the penetration depth are important in 
determining N tum . These local gradients are strongly depen- 
dent upon the details of deleptonization and cooling. There- 
fore, the time-evolution of the characteristic frequencies is 
a strong function of both the dense-matter EOS in the PNS 
below the boundary layer and the local microphysics of the 
boundary layer itself. 

Figure 9 compares f p with the spectrogram, dE^/df 
(eq. 20), of the simulation using the 15-M progenitor and 
L Ve = 3.7 x 10 52 erg s" 1 . From bounce until explosion, the 
frequency at peak power coincides with f p . This strong cor- 
relation persists whether prompt convection, convection, or 
the SASI are the dominant instabilities and implies that the 
plume's buoyant deceleration determines the characteristic 
GW frequency at all times before explosion. 

In Fig. 15, we plot f p as a function of time after bounce for 
all 2D simulations. Each shade (color in the online version) 
represents a set of simulations with the same progenitor model 
but different L Ve : black for 12 M , blue (online) for 15 M , 
green (online) for 20 M©, and red (online) for 40 M@. As in 
Fig. 9, careful inspection of these /y,-time curves show that 
they agree with the trends of peak frequency in the GW spec- 
trograms, Fig. 8. For simulations that use the same progen- 
itor model, their f p evolutions are practically indistinguish- 
able. The only differences are that the higher luminosity sim- 
ulations explode earlier (the end of / p -time curves mark the 
approximate time of explosion) and consequently have lower 
frequencies at explosion. In general, the higher mass progen- 
itors produce higher frequencies, with the exception that the 
12- and 15- M Q models have very similar f,-time curves. 
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Yoshida et al. (2007) also investigated the interaction be- 
tween the plumes of postshock-convection/SASl region and 
the stable layers below it. Their motivation was to investi- 
gate the excitation of PNS core g-modes, which are vital for 
the success of the acoustic mechanism (Burrows et al. 2006, 
2007b), by these plumes. They concluded that the driving 
pressure perturbations had frequencies that are much too low 
to cause excitation of the PNS core g-modes. As a result, they 
suggest that core g-mode excitation is inefficient, rather the 
core g-modes are forced oscillations. However, their analysis 
relied upon the power spectrum of pressure perturbations at 
the base of 2D simulations with inner boundaries at ^50 km. 
While they attempt to handle the subsonic outflow through 
this inner boundary consistently, they are forced to set the ra- 
dial velocity at this inner boundary to maintain the subsonic 
structure above the boundary. This boundary condition cer- 
tainly limits the motions at the base of their calculations, and 
their unfortunate choice of the location of their inner bound- 
ary is exactly where the plume's motions are most important 
for exciting core g-modes. 

Furthermore, there is some confusion in the literature con- 
cerning the appropriate frequencies and excitation mechanism 
for the core g-modes. While pressure perturbations play some 
role in exciting core g-modes, it is well known that buoy- 
ancy is the dominant driving force in g-waves (e.g. Unno et al. 
1989, Ch III, §13). Therefore, velocity perturbations, which 
result from both pressure and buoyancy forces, are more rele- 
vant in characterizing the perturbation spectrum at the surface 
of the PNS. The power spectrum of the pressure perturbations 
is most useful at highlighting the sounds waves that are emit- 
ted by oscillating PNS. Furthermore, in their discussion of the 
relevant frequencies, Yoshida et al. (2007) reference the ~30 
to ~80 Hz shock-oscillation frequencies in Table 1 of Bur- 
rows et al. (2007b) and conclude that motions associated with 
these frequencies are much too low to excite core g-modes 
which typically have frequencies around 300 Hz. However, 
these shock-oscillation frequencies are not relevant in driv- 
ing oscillations at the surface of the PNS. Rather, we have 
shown that the characteristic frequencies associated with the 
plume impulses are most relevant and have frequencies that 
are very similar to the core g-mode frequencies. Though we 
do not address the acoustic mechanism nor the core g-modes 
in this paper, in light of our characteristic frequency analysis, 
it is possible that the plumes associated with postshock con- 
vection/SASI excite core g-mode oscillations (Burrows et al. 
2007b). 

4. DISCUSSION AND SUMMARY 

In this paper, we present a systematic study of the GW 
emission from 2D core-collapse supernova simulations in the 
context of the convection/SASI-aided neutrino mechanism. 
Using a range of progenitor masses (12, 15, 20, and 40 M Q ) 
and neutrino luminosities, we investigate the effects of core 
structure and explosion on the GW signal. In addition to 
GW signatures of prompt convection, PNS and postshock 
convection, and the SASI, we characterize the GW signal of 
asymmetric mass ejections during explosions and identify 
the signatures of explosion that current and near-future GW 
observatories may detect. Moreover, we propose that the 
characteristic GW frequencies and amplitudes are deter- 
mined by the deceleration of postshock-convective/SASI 
plumes by the buoyantly stable layers below. We pro- 
vide gravitational waveforms from all simulated models at 



http : / / stellarcollapse . org/ gwcat a log /murphy 

For all exploding models, we highlight four distinct phases 
in the matter GW strain, h+D. Figure 2 shows that just after 
bounce, prompt convection leads to a burst in GW emission 
with h+D reaching ~5 cm and having a frequency of ~400 
Hz. This lasts for ~50 ms or so before steady-state con- 
vection is established and h+D settles down to ~1 cm. As 
the SASI increases in vigor, h+D increases to ~10 cm and is 
punctuated with strong spikes that are strongly correlated with 
downdrafts striking the PNS "surface." From bounce until ex- 
plosion, the peak frequency of the GW emission rises from 
^100 Hz to ^300-400 Hz, with higher values corresponding 
to higher mass progenitors and/or later explosion times (lower 
neutrino luminosities). 

For a few simulations, the transition to strong SASI and 
GW emission is initiated by a substantial drop in accretion 
rate, which is caused by an abrupt change in density/entropy 
at the base of the progenitor's oxygen-burning shell (Mur- 
phy & Burrows 2008b). The fact that the correlation is weak 
precludes using individual GW detections to investigate the 
location of this entropy edge, but a large sample of detec- 
tions could provide statistical constraints. However, if con- 
vection/SASI is the only GW emission mechanism, detection 
of CCSNe by current and near-future GW detectors is lim- 
ited to our galaxy. Within this radius, the rate of CCSNe is 
roughly one per century (e.g., van den Bergh and Tammann 
1991; Cappellaro et al. 1999), a rate that is far too low to 
build a statistical sample. Such a sample of GW detections 
from CCSNe is only feasible if they can be detected with dis- 
tances out to ^4-5 Mpc, where the rate increases to ~l-2 per 
year (e.g., Ando et al. 2005). Third generation observatories, 
e.g., the proposed Einstein Telescope 13 , or even more sensi- 
tive detectors may be required to study the GW emission from 
convection/SASI in CCSNe out to ~4-5 Mpc. 

As explosion ensues, the SASI and neutrino-driven convec- 
tion above the neutrinosphere subside, causing the rapid os- 
cillations and large spikes in the GW strain to reduce in am- 
plitude. In their stead, aspherical mass ejection results in a 
slow, but sometimes large, monotonic rise in the GW strain. 
For spherical explosions, the strain settles back to zero after 
the onset of explosion, while the strain rises to positive val- 
ues for prolate explosions and drops to negative values for 
oblate explosions. This is similar to the "burst with mem- 
ory" GW signal that asymmetric neutrino emission produces 
(Braginskii & Thorne 1987; Burrows & Hayes 1996; Miiller 
& Janka 1997). Consequently, the peak GW power emitted 
moves from ^300-400 Hz to ~10s of Hz. Initial and En- 
hanced LIGO sensitivity at low frequencies inhibit direct de- 
tection of this explosion signature. Although the improved 
sensitivity of Advanced LIGO at lower frequencies might en- 
able direct detection of the explosion signature, the abrupt re- 
duction in power at the higher frequencies is the most obvious 
signature of explosion for the first generation detectors. 

We stress that this reflects the integrated, £ = 2 morphol- 
ogy of the mass during the explosion. The snapshots of en- 
tropy in Fig. 4 show that while the general morphology of the 
matter may be prolate, oblate, or spherical, the entropy and 
composition distributions are much more complex. So, while 
the shock wave may be spherical, the composition distribu- 
tions are likely to be aspherical. The spherical shock wave 

13 http : / /www . et-gw . eu 
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but asymmetric Fe and Si ejecta distributions in the Cas A su- 
pernova remnant, might be a prime example of this (DeLaney 
et al. 2009). 

These monotonic changes due to asymmetric matter ejec- 
tion will be superposed on the large secular changes due 
to asymmetric neutrino emission (Burrows & Hayes 1996; 
Miiller & Janka 1997; Miiller et al. 2004; Marek et al. 2009). 
If the timescales and amplitudes are similar, then it will be 
difficult to distinguish between the neutrino and asymmetric 
matter ejection GW signals. Simulations that employ more 
consistent neutrino transport for several 100 ms past explosion 
are necessary to resolve this issue. Regardless, for CCSNe 
in which convection/SASI is the dominant source of GWs, 
distinguishing among prolate, oblate, or spherical explosions 
will require direct detection of the GW amplitude to enable 
a determination of the sign of the amplitude. This will only 
be possible if a very close post-main-sequence star such as 
Betelgeuse, which is a mere 197 ±45 pc away (Harper et al. 
2008), explodes in the near future. Otherwise, the distances 
to typical Galactic SNe (^10 kpc), detector characteristics, 
and their sensitivity to frequencies >50 Hz imply that the low 
frequency "memory" signal will not be detected, and the sig- 
nature of explosion will be the abrupt reduction in GW power. 

Except for the explosion signature, Marek et al. (2009) re- 
port similar characteristics in the GW strain. The evolution 
in frequency from ^100 Hz just after bounce to larger fre- 
quencies at later times agrees with the results of Marek et al. 
(2009). In fact, our higher frequencies are consistent with the 
calculations that use a stiff EOS, but roughly half those re- 
ported for their softer EOS simulation. Kotake et al. (2007) 
and Kotake et al. (2009), on the other hand, report only one 
phase of weak GW emission exhibiting gradual growth and 
one characteristic frequency at ^100 Hz. We suspect that the 
idealized initial conditions and omission of the inner 50 km 
precludes them from simulating many of the features that we 
and Marek et al. (2009) identify. 

While the total GW emission originates from regions in- 
cluding the PNS convection at depth, the stable layer above 
it, and the base of the outer convective and SASI regions, 
we find that the characteristic GW frequencies and ampli- 
tudes are, for the most part, established by the interaction of 
postshock-convection/SASI region and the stable layers be- 
low. As plumes in the convective/SASI region fall toward the 
stable layer below, momentum carries their motion beyond the 
lower boundary of convection as defined by the Ledoux con- 
dition. As they penetrate the stable regions, the plumes expe- 
rience a buoyancy acceleration upward. It is this impulse that 
gives rise to the GW emission. The pulse width has a charac- 
teristic frequency, f p , that corresponds to the peak power of 
GW emission, and we show that f p is insensitive to the plume 
velocity. Instead, f p is set by the buoyancy frequency where 
the plume reverses its downward motion, f p ~ AWn/(27r). 
Moreover, the conditions during this deceleration are of the 
correct order of magnitude to explain the GW strain ampli- 
tude. In particular, we show that the amplitude is directly 
proportional to the plume velocity, h+ oc f p v p oc N tum v p . 

Since f p is set by Mum, the frequency at peak GW power 
gives direct information on the core structure. This is not too 
surprising. In fact, using soft and stiff dense-matter EOSs, 
Marek et al. (2009) find higher peak frequencies for the for- 
mer, which they attribute to the compactness of the PNS. 
While compactness certainly influences N tum , we emphasize 



that local conditions at the PNS "surface" and their evolution 
are just as important in determining the characteristic GW fre- 
quencies. As a point of clarification, we find that sloshing 
and turnover timescales do not set the characteristic GW fre- 
quency. Assuming that / /; is proportional to v p divided by a 
length scale such as an eddy size, one might conclude that f p 
is directly proportional to the plume speed. This contradicts 
our earlier finding that f p is dependent only on the core struc- 
ture. However, in §3.4, we show that the relevant length scale 
is the deceleration length, or penetration depth, D p , and that 
D p ~ v p /N, giving the result that f p is proportional to N and 
independent of v p . 

In general, we find that more massive progenitors have 
higher f p at a given time, though the 12- and 15-M Q mod- 
els show similar f p evolution. Interestingly, our simulations 
show that the f p vs. time evolution for a given progenitor 
model is independent of the neutrino luminosity. The fact that 
the characteristic frequency at explosion is dependent upon 
neutrino luminosity is simply due to the fact that higher lu- 
minosity simulations explode earlier in an otherwise similar 
/ ; ,-time evolution. 

Though we state that our conclusions are qualitatively valid, 
we have made several assumptions that enabled timely calcu- 
lation of the simulations in this paper. The effects of neu- 
trino transport were approximated by local heating and cool- 
ing functions, azimuthal symmetry was assumed, though na- 
ture is inherently 3D, and we used spherical Newtonian grav- 
ity rather than a full general relativistic treatment. 

Differences between 2D and 3D simulations are likely to 
matter quantitatively. For one, in steady-state stellar convec- 
tion, the eddies are larger and have faster speeds in 2D sim- 
ulations (Arnett et al. 2009; Meakin & Arnett 2007b, and the 
references therein). Although we have demonstrated that the 
characteristic frequencies are unlikely to be affected, the GW 
strain amplitude is directly proportional to the plume speed 
and will be affected. However, these arguments are for con- 
vection alone; they ignore the fact that in the gain region, 
convection, the SASI, and advection are equally important. 
For example, the fact that we can explain the characteristic 
GW frequencies by considering the interaction of convective 
plumes with neighboring buoyant layers implies that convec- 
tion is important. On the other hand, the rapid rise in plume 
velocities as the SASI becomes nonlinear is a clear indica- 
tion that the SASI is equally important. Even though simu- 
lations are beginning to explore the 3D aspects of the SASI 
(Iwakami et al. 2008b, a; Kotake et al. 2009) and advection- 
modified convection has been explored in the linear regime 
(Foglizzo et al. 2006), there has yet to be a complete analy- 
sis of nonlinear convection that is modified by either. Fur- 
thermore, Murphy & Burrows (2008b) have noted that the 
accretion through the convective/gain region imposes an in- 
herent asymmetry in the convective flow. While some matter 
advects quickly through the convective/gain region (these are 
the downdrafts/plumes), some lingers for very long times. To 
what degree do these effects hold in 3D simulations? As yet 
the similarities and differences are not obvious and remain an 
interesting and important avenue of research for core-collapse 
simulations. 

Though 3D effects may not alter the characteristic frequen- 
cies, including general relativity and consistent neutrino trans- 
port undoubtedly will. With general relativistic corrections, 
the PNS structure is more compact, leading to higher char- 
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acteristic frequencies. Without proper treatment of neutrino 
transport, in particular, and v r emission, the PNSs in our 
simulations do not lose energy and shrink as they should. In- 
cluding these effects should increase (GM,/^) 1 ' 2 over time, 
further increasing f p in our simulations. Furthermore, the 
neutrino transport is important in determining the structure, 
and, hence, buoyancy frequency, at the surface of the PNS. 
Therefore, to obtain a reliable and quantitative mapping be- 
tween the model characteristic frequencies from postshock- 
convection/SASI motions and those observed by GW detec- 
tors such as LIGO, approximations to GR and the neutrino 
transport that more faithfully mimic nature are required in fu- 
ture simulations. 
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FIG. 1 . — The dimensionless GW strain h+ times the distance, D, vs. time after bounce for the suite of simulations presented in this paper. Each panel represents 
a single progenitor model and the curves are labeled by L Ve in units of 10 52 erg s . The accretion rates (in Mq s -1 ) at the stalled shock of non-exploding models 
are shown for comparison. All models show GW features associated with prompt convection, postshock convection, SASI, and explosion with the exception that 
the run using M = 40 Mq and L Vc = 6.0 X 10 52 erg s _1 does not explode. See Fig. 2 for a sample GW strain evolution with these features labeled. In general, 
lower luminosity runs develop the strong SASI motions and explosions at later times, but all runs saturate at roughly similar amplitudes (h + D ~ 10 cm). In the 
simulations using the 2O-M0 progenitor, the accretion of the interface between the progenitor's Si and O burning shells through the shock causes a sudden drop 
in accretion rate. This initiates strong convective/SASI motions in the highest luminosity runs (3.4, 3.6, and 3.8) and suggests a means to use GW signals to probe 
the location of this interface. However, this feature is unique to the runs using the 2O-M0 progenitor. 
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FIG. 2. — A sample of GW strain (h+) times the distance, D, vs. time after bounce. This signal was extracted from a simulation using a 15-Mq progenitor 
model (Woosley & Heger 2007) and an electron-type neutrino luminosity of L Ve = 3.7 X 10 52 erg s~'. Prompt convection, which results from a negative entropy 
gradient left by the stalling shock, is the first distinctive feature in the GW signal from to ~50 ms after bounce. From ~50 ms to ~550 ms past bounce, the 
signal is dominated by PNS and postshock convection. Afterward and until the onset of explosion (~800 ms), strong nonlinear SASI motions dominate the 
signal. The most distinctive features are spikes that correlate with dense and narrow down-flowing plumes striking the "PNS" surface (~50 km). Around ~800 
ms the model starts to explode. In this simulation, the GW signal during explosion is marked by a significant decrease in nonlinear SASI characteristics. The 
aspherical (predominantly prolate) explosion manifests in a monotonic rise in h+D that is similar to the "memory" signature of asymmetric neutrino emission. 
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FIG. 3. — Plots of h+D vs. time after bounce for three models showing the effect of global asymmetries in explosion. The general shapes of the explosions are 
prolate (M = 12M Q ,L, /( , = 3.2 X 10 52 erg s" 1 , and thick light-gray line), spherical (M = 12 Mq, L Ve =2.2 X 10 52 erg s~\ and thin black line), and oblate (M = 20 
MQ,i„ ( = 3.4 X 10 52 erg s , and thin gray line), and correspond to positive, zero, and negative "memory," respectively. See Fig. 4 for color maps of the entropy 
distribution during explosion for the same models. 
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FIG. 4. — Snapshots of the entropy (Ag/baryon) distribution during explosion for the three models highlighted in Fig. 3. Lighter shades (wanner colors in the 
online version) represent higher entropies, while darker shades (cooler colors in the online version) shades represent lower entropies. The global asymmetry of 
the matter interior to the shocks is oblate (M = 20 Mq and L Ve = 3.4 X 10 52 erg s~'), prolate (M = 12 M Q and h Ve = 3.2 X 10 52 erg s" 1 ), and spherical (M = 12 
Mq.L^ = 2.2 X 10 52 erg s~'). Even though the GW signal may indicate that the explosion is in general "spherical," "oblate," or "spherical," the entropy (hence 
temperature and composition) distributions may show higher order structure. 
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FIG. 5. — Gray-scale (color in the online version) map of the integrand that leads to the GW signal (See eqs. 12, 13, and 21). Lighter shades (brighter colors in 
the online version) correspond to higher values. The results are for the 15-Mq progenitor and L Ve = 3.7 X 10 52 erg s . The left panel represents the stalled-shock 
phase during vigorous postshock convection/SASI motions. During this phase, the signal originates predominantly from PNS convection and the deceleration of 
plumes below the gain radius (~100 km). The right panel shows that the "memory" signal during explosion is a result of acceleration at an asymmetric shock. 
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FIG. 6. — GW waveforms, h+D vs. time, showing the contributions of PNS convection and the SASI. The model shown is the same as in Fig. 2. For reference, 
the entire GW signal is shown in both panels (solid-black line). The signal originating from > 50 km (orange, online version) and < 50 km (blue, online version) 
is shown in the top panel. This radius is roughly the division between nonlinear SASI motions and PNS convection motions. Most, but not all, of the signal 
associated with prompt convection originates in the outer convection zone. There is a non-negligible contribution from the PNS. The monotonic rise in the GW 
strain during explosion clearly originates from the outer (exploding) regions. Even though the region for the convection/SASI and its nonlinear motions are above 
50 km, these motions influence the PNS convective motions below 50 km. It is telling that once the model explodes, and the nonlinear SASI motions subside but 
the PNS convection does not, the GW signal from below 50 km diminishes as well. The bottom panel shows the GW signal from five regions, each with different 
outer radii (30, 40, 50, 60, and 100 km). The strengthening of the GW signal associated with the SASI is apparent for all, suggesting that the influence of the 
SASI diminishes only gradually with depth. 
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FIG. 7. — /! c har (eq. 17) vs. frequency for the suite of simulations presented in this paper. The spectra show broad peaks and some dependence upon the 
progenitor mass: ~300 Hz for 12 Mq and ^400 Hz for 40 Mq. The power of the spectra show only a slight increase (if at all) with neutrino luminosity. For 
comparison, the approximate noise thresholds for Initial LIGO (solid-black curve), enhanced LIGO (dot-dashed-black curve) and Advanced LIGO (dashed-black 
curve) are plotted. 
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TABLE 1 

Quantitative Summary of GW Characteristics' 



Mass^ 




\h I 4 
'fl+.max | 


S/Ni irn 5 


S/N.i irr, 6 


Lj ' 1 'aavLKjU 


h,.w„ * 

"cnai ,max 


J cnar,max 




12 


1.8 


0.27 


0.8 


2.0 


7.9 


1.04 


253 


6.1 


12 


2.2 


0.20 


0.7 


1.7 


6.7 


0.61 


200 


3.5 


12 


2.8 


0.46 


1.2 


2.7 


10.8 


1.35 


273 


8.5 


12 


3.2 


1.13 


1.6 


4.0 


15.7 


2.23 


294 


20.4 


15 


3.2 


0.37 


1.1 


2.6 


11.4 


2.14 


345 


17.5 


15 


3.4 


0.40 


1.0 


2.4 


10.2 


1.49 


406 


14.4 


15 


3.7 


0.37 


1.1 


2.8 


11.5 


1.83 


365 


12.8 


15 


4.0 


1.53 


2.1 


5.3 


21.2 


3.10 


347 


46.1 


20 


3.2 


0.32 


1.4 


3.4 


14.5 


3.00 


348 


29.5 


20 


3.4 


0.70 


1.8 


4.5 


19.5 


4.68 


347 


57.8 


20 


3.6 


0.56 


1.5 


3.9 


16.1 


2.67 


429 


34.5 


20 


3.8 


0.48 


1.5 


3.8 


15.7 


3.14 


369 


33.8 


40 


6.0 


0.33 


1.0 


2.6 


12.2 


3.23 


423 


36.2 


40 


10.0 


0.68 


1.4 


3.6 


17.3 


3.93 


359 


47.1 


40 


12.0 


0.82 


1.4 


3.4 


14.0 


2.04 


323 


16.6 


40 


13.0 


4.53 


0.8 


1.8 


9.4 


1.08 


It 


4.1 



This table lists the integrated GW characteristics of the 2D simulations. These simulations represent a two-dimensional 
parametrization that investigates the dependence of GW emission on progenitor mass (column 1) and neutrino luminosity 
(column 2). 

^Progenitor model (Mq). 

"^Neutrino Luminosity (10 52 erg s -1 ). 

4 Maximum GW strain (1(T 21 at 10 kpc). 

^Optimal theoretical single-detector signal-to-noise using the initial LIGO sensitivity curve (Gustafson et al. 1999). 
^Optimal theoretical signal-to-noise using the Enhanced LIGO sensitivity curve (Adhikari 2009). 
^Optimal theoretical single-detector signal-to-noise using the burst-mode Advanced LIGO sensitivity curve (Shoemaker 
2006). 

^Maximum of the characteristic strain spectrum defined in eq. 17 (10~ 21 at 10 kpc). 
^Frequency location of /i € har.max (Hz). 
10 GW energy emitted (10" n M c 2 ). 

^Because this simulation explodes early and with large asymmetry, the low frequency "memory" signature in the GW 
strain dominates the energy spectrum. 
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FIG. 8. — Spectrograms of the GW signals showing dEX^/df as a function of frequency (vertical axes) and time (horizontal axes) after bounce for the entire set 
of simulations presented in this paper. Bursts of power are associated with prompt convection and the SASI. Frequency at peak power increases from ~ 100 Hz to 
^300-400 Hz, depending upon the explosion time and progenitor used. We show in §3.4 that this frequency and the trend to higher frequencies is a consequence 
of the core structure. During explosion, GW power transitions to lower frequencies (~10s of Hz). 
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Time ofter bounce [s] 

FIG. 9. — Spectrogram, dE^/df, and characteristic plume deceleration frequency, /,,, vs. time after bounce for the simulation using the 15-Mq progenitor and 
Li/ e =3.1 X 10 52 erg s . Similar to Fig. 2, the three features showing the largest power are prompt convection just after bounce, nonlinear SASI plumes/motions, 
and explosion. Unlike in the GW strain, the specific morphology of the explosion, other than it is asymmetric, is not discernible in the spectrogram. From bounce 
until explosion, the frequency at peak power agrees with f p . This strong correlation persists whether prompt convection, convection, or the SASI are the dominant 
hydrodynamic processes and implies that the plume's buoyant deceleration determines the characteristic GW frequency. 
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FIG. 10. — Comparison of h+D with the maximum downward plume speed below 120 km, v p , vs. time after bounce. The increase in GW wave strain at ^600 
ms coincides with the rise in v p , and many spikes in h+D coincide with rapid changes in v p . A rough estimate of h+D using these plume speeds gives ~5 cm, 
which is consistent with the calculated GW strain (see the text in §3.4). 
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FIG. 11. — The buoyancy frequency, N (blue-solid line, online version), (GM, /r 3 ) 1 / 2 (purple-dashed line, online version), and the radial velocity, v, (black 
points), vs. radius at 300 ms after bounce for the 2D simulation that uses the 15-Mq progenitor model and L Ve = 3.7 X 10 52 erg s _1 . Starting at larger radii, 
the velocity points show spherical infall, an asymmetric shock, postshock convection/SASI, a stable layer, the PNS convection, and a stable inner core. At this 
time the dynamics are dominated by postshock convection and mild SASI oscillations. The convection regions are characterized by two regions: 1) that which is 
unstable by the Ledoux condition (green lines, online version) and 2) the region of overshoot below the postshock-convection region (orange line, online version). 
The penetration depth, D p , associated with overshoot is computed by integrating the plume's equation of motion and finding the depth at which the plume turns 
around. The plume frequency, f p , obtained from the HWHM timescale of the acceleration pulses, is ~ 196 Hz. 
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FIG. 12. — Similar to Fig. 11, except that t = 620 ms. At this time, the SASI oscillations are much stronger, resulting in a higher penetration depth of D p ~ 39 
km. Because N has also increased, the characteristic frequency, f p , is higher at 272 Hz. 
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FIG. 13. — Shown are characteristic plume frequency, /,,, the buoyancy frequency at Dp, AWn, and (GM,-/P) [ / 2 at Dp as a function of time after bounce for 
the 2D simulation using the 15-Mq progenitor and L Ue = 3.7 X 10 52 erg s . Initially, all three coincide. However, while AWn tracks the time-evolution of f p , 
(GMr/r 3 ) 1 / 2 diverges from the other two measures. This indicates that while (CM r /r 3 )'' 2 (which is related to the compactness of the PNS) sets the general 
frequency scale, the local logarithmic gradients at the penetration depth are just as important in determining AWn- 
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FIG. 14. — The ratio N tmn D p /v p , which is the square root of Rj, after inserting the approximation that Ab ~ N 2 D P , vs. time after bounce for all 2D simulations 
listed in Table 1 . Each line is labeled by the progenitor mass and neutrino luminosity used. Each shade (color in the online version) represents a set of simulations 
using the same progenitor model: black for 12 Mq, blue (online) for 15 Mr, green (online) for 20 Mq, and red (online) for 40 Mq. The style of line corresponds 
to the neutrino luminosity used, but note that the range of neutrino luminosities for a set of models that use the same progenitor mass is different from another 
set. This ratios are ~ 1.5 to ~2, and change less than 30% during the simulations. This validates both the assumption that Rt is of order 1 where the plumes turn 
around and the approximation that Ab ~ N 2 D P . 
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FIG. 15. — Characteristic plume frequency, as a function of time after bounce for all 2D simulations listed in Table 1. The shading (color in the online 
version) and line schemes are the same as in Fig. 14. For simulations that use the same progenitor model, the / p -time evolutions are close. The only differences 
are that the higher luminosity simulations explode earlier (the end of / p -time curves mark the approximate time of explosion) and consequently have lower 
frequencies at explosion. In general, the higher mass progenitors produce higher frequencies, with the exception that the 12- and 15-Mq models have very 
similar / p -time curves. 



